Global Optimization by Basin-Hopping and the Lowest Energy Structures of 
Lennard- Jones Clusters Containing up to 110 Atoms 
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We describe a global optimization technique using 'basin-hopping' in which the potential energy 
surface is transformed into a collection of interpenetrating staircases. This method has been designed 
to exploit the features which recent work suggests must be present in an energy landscape for efficient 
relaxation to the global minimum. The transformation associates any point in configuration space 
with the local minimum obtained by a geometry optimization started from that point, effectively 
removing transition state regions from the problem. However, unlike other methods based upon 
hypersurface deformation, this transformation does not change the global minimum. The lowest 
known structures are located for all Lennard- Jones clusters up to 110 atoms, including a number 
that have never been found before in unbiased searches. 
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I. INTRODUCTION 

al optimization is a subject of intense current 
interest .La Improved global optimization methods could 
be of great economic importance, since improved solu- 
tions to travelling salesman type problems, the routing of 
circuitry in a chip, the active structure of a biomolecule, 
etc., equate to reduced costs or improved performance. In 
chemical physics the interest in efficient global optimiza- 
tion methods stems from the common problem of finding 
the lowest energy configuration of a (macro)molccular 
system. For example, it seems likely that the native 
structure of a protein is structurally related to the global 
minimum of its potential energy surface (PES). If this 
global minimum could be found reliably from the pri- 
mary amino acid sequence, this knowledge would pro- 
vide new insight into the nature of protein folding and 
save biochemists many hours in the laboratory. Unfortu- 
nately, this goal is far from being realized. Instead the 
development of global optimization methods has usually 
concentrated on much simpler systems. 

Lennard- Jones (LJ) clusters represent one such test 
system. Here the potential is 
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where e and 2 1 / 6 <r are the pair equilibrium well depth and 
separation, respectively. We will employ reduced units, 
i.e. e = a = 1 throughout. Much of the initial inter- 
est in LJ clusters was motivated by a desire to calculate 
nucleation rates for noble gases. However, as a result 
of the wealth of data generated, the LJ potential has 
been used not only for studying global optimization but 
also the effects of finite size on phase transitions such as 
melting. Through the combined efforts of many workers 
likely candidates for the global minima, of LJjy clusters 
have been found up to N — 147.oEil This represents 



a significant achievement since extrapolation of Tsai and 
Jordan's comprehensive enumeration of minima for small 
LJ clusterstH suggests that the PES of the-147-atom clus- 
ter possesses of the order of 10 60 minima.E3 

Previous -studies have revealed that the Mackay 
icosahedrorJlj provides the dominant structural motif for 
LJ clusters in the size range 10-150 atoms. Complete 
icosahedra are possible at N = 13, 55, 147, ... At most in- 
termediate sizes the global minimum consists of a Mackay 
icosahedron at the core covered by a low energy overlayer. 
As a consequence of the phase behaviour of LJ clusters, 
finding these global minima is relatively easy. Studies 
have shown that in the region of the solid-liquid transi- 
tion the cluster is observed to change back and forth_bj£=. 
tween a liquid-like form and icosahedral structuresBj'EJ 
As a result of this 'dynamic coexistence' a method as 
crude as molecular dynamics within the melting region 
coupled with systematic minimization of configurations 
generated by the trajectory is often sufficient to locate 
the global minimum E3 




FIG. 1. Non-icosahedral Lennard- Jones global minima. 
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However, there are a number of sizes at which the 
global minimum is not based on an icosahedral struc- 
ture. These clusters are illustrated in Figure For LJ 38 
the lowest energy structure is a face-centred-cubic (fee) 
truncated octahedron,BM and for N = 75, 76, 77_J02, 
103 and 104 geometries based on Marks' decahedral are 
lowest in energy.t-Ml3 For these cases finding the lowest 
minimum is much harder because the global minimum 
of free energy only becomes associated with the global 
potential energy minimum at temperatures well below 
melting where the dynamics of structural relaxation are 
very slow. For LJ38, the microcanonical temperature 
for the transition from face-centred cubic to icosahe- 
dral structures has been estimated to be about 0.12efc _1 , 
where k is the Boltzmann constant, and for LJ75 the esti- 
mate for the decahedral to icosahedral transition is about 
QMek^H (For comparison, melting typically occurs at 
about T = 0.2- 0.3efc _1 .) 

The topography of the PES can also play a kgjf role 
in determining the ease of global optimizations. A 
detailed study of the LJ 38 PES has shown that there 
is a large enecgy barrier between the fee and icosahe- 
dral structuresc3 which correspond to well-separated re- 
gions of the PES. Furthermore, fee and decahedral struc- 
tures have less polytetrahedral character than icosahe- 
dral structures, and hence they have less in common with 
the liquid-like state, which. is-jcharacterized by disordered 
polytetrahedral packing.EZrEj Since the vast majority of 
configuration space is dominated by 'liquid-like' configu- 
rations, it is therefore harder to find global minima based 
upon fee and decahedral packing using unbiased searches. 

These considerations explain why global optimization 
methods have only recently begun to find the truncated 
octahedroriOESEJO and why, until now, the Marks' dec- 
ahedron has never been found by an unbiased global op- 
timization method. The greater difficulty of finding the 
LJ75 global minimum compared to LJ38 can probably be 
explained by the slightly smaller transition temperature, 
the sharper transition and the much larger number of 
minima on the LJ75 PES. 

Before we consider the effectiveness of different global 
optimization methods for Lennard-Jones clusters it is 
interesting to note that the use of physical principles 
to consiructp-eiiod candidate structures for the global 
minimaoTjElO'Ej or to reduce the configuration space 
that needs to be searchedEJHLd led to the initial discov- 
ery of 93% of the LJ global minima in this size range. It 
seems that physical insight into a specific problem will of- 
ten be able to beat unbiased global optimization, a view 
expressed by Ngo et alS3 in their discussion of computa- 
tional complexity. 

One difficulty in evaluating the relative performances 
of different global optimization methods is that, too of- 
ten, the methods have only been applied to small clus- 
ters, or to larger clusters with global minima that are 
especially stable, such as LJ55. It is also difficult to 
draw any firm conclusions about how efficient different 
methods may be when the number of searches employed 



varies widely. However, it seems reasonable to suggest 
two hurdles that any putative global optimization ap- 
proach should aspire to. The first is the location of the 
truncated octahedron for LJ38, and any method which 
fails this test is unlikely to be useful. The second hur- 
dle is the location of the Marks decahedron for LJ75; 
this problem poses a much more severe test for an unbi- 
ased search and one which does not appear to have been 
passed until the present work. 

The most successful global optimization results 
for LJ clpjtei»-| reported to date are for genetic 
algorithms.li§l^H'E-3 These methods mimic some aspects 
of biological evolution: a population of clusters evolves 
to low energy by mutation and mating of structures, and 
selection of those with low potential energy. To be suc- 
cessful new configurations produced by 'genetic manip- 
ulation' are mapped onto minima by a local optimiza- 
tion algorithm such as the conjugate gradient method. 
The study by Deaven et al. is particularly impressive, 
since these workers matched or beat all the lowest en- 
ergy minima that they knew of up to N — 100, includ- 
ing the truncated octahedron (although they probably 
missed the global minima at A^=69 and 75-78). Niesse 
and Mayne were also able to locate the LJ38 truncated 
octahedron, and report that this structure took about 25 
times longer to find than the icosahedral global minima 
of the neighbouring sizes. 

Another class of global optimization techniques, some- 
times called hypersurface deformation methods, at- 
tempts to simplify the problem by applying a transfor- 
mation to the PES whiph^sinoothes it and reduces the 
number of local minimaO'Lj The global minimum of the 
deformed PES is then mapped back to the original sur- 
face in the hope that this will lead back to the global min- 
imum of the original PES. The distinctions between the 
various methods of this type lie in the type of transforma- 
tions that-ape used, which include applying the diffusion, 
equationj23l2a increasing the range of the potentialj23'c2l 
and shifting thepppsition of the potential minimum to- 
wards the originH The performance of hypersurface de- 
formation methods has been variable: the potential shift 
approach managed to find the 38-atom truncated octa- 
hedron, but other workers report difficultiesEJEj for the 
trivial cases of LJg and LJ9 where there are only 8 and 
21 minima on the PES, respectively. 

Although intuitively appealing, the problem with hy- 
persurface deformation is that there is no guarantee that 
the global minimum on the deformed PES will map onto 
the global minimum of the original surface. This dif- 
ficulty is clearly illustrated when we consider Stillingcr 
and Stillinger's suggestion of increasing the range of the 
potential:!^] it has been shown that the global minimum 
may in fact depend rather sensitively on the range of 
the potential, with t he a ppearance of numerous 'range- 
induced' transitions.t3c3 

Other methods include those based on 'annealing'. 
Such approaches take advantage of the simplification in 
the free energy landscape that occurs at high tempera- 
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tures, and attempt to follow the free energy global mini- 
mum as the temperature is decreased. At zero Kelvin the 
free energy global minimum and the global minimum pf, 
the PES must coincide. Standard simulated annealingo 
was .used by Wille to find a few new minima at small 
sizesQ but does not appear to have been systematically 
applied to LJ clusters. More sophisticated variants of 
this technique-, include gaussian density annealing and 
analogues JijacH but again some appear to fail at small 
sizes O'cZl 

The difficulty with the annealing approach methods 
is that, if the free energy global minimum changes at 
low temperatures where dynamical relaxation is slow, 
the algorithms will become stuck in the structure cor- 
responding to the high temperature free energy global 
minimum. Such methods are therefore likely to experi- 
ence difficulties in finding the global minima for LJ38 and 
LJ75. In the language employed in recent protein folding 
literature p3 annealing will fail when Tf < T g , where Tf 
is the 'folding' temperature below which the global po- 
tential energy and free energy minima coincide, and T g 
is the 'glass' temperature at which the system effectively 
becomes trapped in a local minimum. 

Another method which attempts to reduce the effects 
of barriers on the PES makes use of quantum tunnelling. 
The diffusion Monte Carlo approach is used to find the 
ground state wavefunction, which should become locsA-, 
ized at the global minimum as h is decreased to zero.Ea 
A more rigorous approach has been applied by Maranas 
and Floudas, who found upper and lower bounds for 
the energy of the global minimum. However, the com- 
putational expense of this method, which scales as 2 N 
with the number of atop^-means that it has only been 
used for small systems .Ed'ELlI Most of the above studios, 
along with the recently described 'pivot method'EHTE-3 
and 'taboo search' J£j't3 have yet to prove their useful- 
ness by passing the first hurdle for LJ 38 suggested above. 
However, this does not necessarily mean that these ap- 
proaches should be discounted, since some authors have 
only applied their algorithms to smaller clusters and may 
not have run enough searches to achieve convergence. 

In the present work we present the results of a 'basin- 
hopping' global optimization technique for Lcnnard- 
Jones clusters. All the known lowest energy structures 
up to N = 110 have been located successfully, includ- 
ing three minima not previously reported. (See Tables 
1 and 2.) This method is also the first unbiased algo- 
rithm to find the global minima based on the Marks 
decahedron around LJ75 and LJio2- For reference, we 
collect the rather scattered results previously reported 
for LJ clusters to provide a complete catalogue of the 
energies and point groups of the lowest energy minima 
that we know of. The results have been collected in 
the first entry of the Ca mbridge Cluster Database at 
http: / /brian. ch.cam.ac.uk . 



II. METHOD 

The present approach has been guided by previous 
work on energy landscapes which has identified features 
that enable the system to locate its global minimum 
efficiently^ In particular, analysis of model energy land- 
scapes using a master equation approach for the dy- 
namics, has provided good evidence that such a surface 
should have a large potential energy gradient and the 
lowest possible transition state energies or rearrangement 
barriers E3 These results immediately suggest a simple 
way to transform the PES which does not change the 
global minimum, nor the relative energies of any local 
minima. We consider the transformed energy E defined 
by: 

E(K)=min{E(K)}, (2) 

where X represents the 3iV-dimensional vector of nuclear 
coordinates and min signifies that an energy minimiza- 
tion is performed starting from X. In the present work 
energy minimizations were performed using the PolaJfc. 
Ribiere variant of the conjugate gradient algorithmic 
Hence the energy at any point in configuration space 
is assigned to that of the local minimum obtained by 
the given geometry optimization technique, and the PES 
is mapped onto a set of interpenetrating staircases with 
plateaus corresponding to the set of configurations which 
lead to a given minimum after optimization. A schematic 
view of the staircase topography that results from this 
transformation is given in Figure ||. These plateaus, or 
basins of attraction, have been visualized in previous 
work as a means to compare the efficiency of different 
transition state searching techniques P 8 " 59 ! 
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FIG. 2. A schematic diagram illustrating the effects of 
our energy transformation for a one-dimensional example. 
The solid line is the energy of the original surface and the 
dashed line is the transformed energy, E. 

The energy landscape for the function E(K) was ex- 
plored using a canonical Monte Carlo simulation at a 
constant reduced temperature of 0.8. At each step all 
coordinates were displaced by a random number in the 
range [—1,1] times the step size, which was adjusted to 
give an acceptance ratio of 0.5. The nature of the trans- 
formed surface allowed relatively large step sizes of be- 
tween 0.36-0.40. For each cluster in the range considered 
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seven separate runs were conducted. Five of these each 
consisted of 5000 Monte Carlo steps starting from differ- 
ent randomly generated configurations of atoms confined 
to a sphere of radius 5.5 reduced units. The subsequent 
geometry optimizations employed a container of radius 
one plus the value required to contain the same volume 
per atom as the fee primitive cell. The container should 
have little effect on any of our results and is only required 
to prevent dissociation during the conjugate gradient op- 
timizations. 

The convergence criterion employed for the conjugate 
gradient optimizations used in the Monte Carlo moves 
need not be very tight. In the present work we required 
the root-mean-square (RMS) gradient to fall below 0.01 
in reduced units and the energy to change by less than 
0.1 e between consecutive steps in the conjugate gradient 
search. Initially it appeared that a tolerance of 0.1 for 
the RMS gradient was satisfactory, but this was subse- 
quently found to cause problems for clusters containing 
more than about 60 atoms. The lowest energy structures 
obtained during the canonical simulation were saved and 
reoptimized with tolerances of 10~ 4 and 10~ 9 for the 
RMS force and the energy difference, respectively. The 
final energies are accurate to about six decimal places. 

Several other techniques were employed in these cal- 
culations, namely seeding, freezing and angular moves. 
Here we used the pair energy per atom, E{i) 1 defined as 
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so that the total energy is 

E=lj2 E (^ (4) 



were taken directly between minima using eigenvector- 
following to calculate pathways .E3 The latter method 
is similar_tp that described recently by Barkema and 
MousseauEj in their search for well-relaxed configurations 
in glasses. Although the computational expense of tran- 
sition state searches probably makes this method uncom- 
petitive for global optimization, our study illustrated the 
possible advantage of working in a space in which only 
the minima are present. The basin- hopping algorithm 
differs in that it is applied in configuration space to a 
transformed surface, rather than in a discrete space of 
minima, and steps are taken stochastically. The genetic 
algorithms described by Deaven et altB and Niesse and 
MayneE9 used conjugate gradient minimization to refine 
the local minima which comprise the population of struc- 
tures that are evolved in their procedure. Hence these au- 
thors are in effect studying the same transformed surface 
as described above, but explore it in a rather different 
manner. We suspect that the success of their methods is 
at least partly due to the implicit use of the transformed 
surface E. 

The transformation of the PES also reduces the bar- 
riers to dissociation. Therefore, to prevent evaporation 
either the cluster can be by placed in a tight-fitting con- 
tainer or the coordinates of the current point in config- 
uration space can be reset to that of the minimum after 
each successful step. In this paper we used the latter 
method and then the present approach is essentially the 
same as the 'Monte Carlo-minimization' algorithm of Li 
and Scheragaj£il who applied it to search the conforma- 
tional space of the pentapeptide [Met 5 ]enkephalin. A 
similar method has recently been used by Baysal and 
Meirovitch£3 to search the conformational space of cyclic 
polypeptides. 

69 78 107 



If the highest pair energy rose above a fraction a of the 
lowest pair energy then an angular move was employed 
for the atom in question with all other atoms fixed, a was 
adjusted to give an acceptance ratio for angular moves 
of 0.5 and generally converged to between 0.40 and 0.44. 
Each angular displacement consisted of choosing random 
9 and 4> spherical polar coordinates for the atom in ques- 
tion, taking the origin at the centre of mass and replacing 
the radius with the maximum value in the cluster. 

The two remaining runs for each size consisted of only 
200 Monte Carlo steps starting from the global minima 
obtained for the clusters containing one more and one less 
atom. When starting from LJ the N — 1 atoms were 
frozen for the first 100 steps, during which only angular 
moves were attempted for the remaining atom, starting 
from a random position outside the core. When starting 
from LJtv +1 the atom with the highest pair energy, E(i), 
was removed and 200 unrestricted Monte Carlo moves 
were attempted from the resulting geometry. 

The above basin-hopping algorithm shares a common 
philosophy with our previous approach in which steps 






FIG. 3. Lennard- Jones global minima that have not pre- 
viously been reported. 



III. RESULTS 

The basin-hopping algorithm has successfully located 
all the lowest known minima up to LJuo, including all the 
non-icosahedral structures illustrated in Figure @ (sizes 
38, 75, 76, 77, 102, 103 and 104) and three new geome- 
tries based upon icosahedra illustrated in Figure || (sizes 
69, 78 and 107). We believe that this is the first time any 
of the six decahedral global minima have been located 
by an unbiased algorithm. The total number of searches 
was fixed in our calculations to provide a simple refer- 
ence criterion. In fact, most of the global minima were 
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found in more than one of the separate Monte Carlo runs. 
The global minima for the smallest clusters were located 
within a few steps in each of the seven runs. To give a 
better idea of how the algorithm performed we will pro- 
vide some more details for the sizes with non-icosahedral 
or newly discovered icosahedral global minima. 

For LJ38 the truncated octahedron was found in four 
out of five of the longer unseeded runs; the first suc- 
cess occurred within a thousand Monte Carlo steps on 
average. Not surprisingly, the global minimum was not 
located in the shorter runs starting from the structurally 
unrelated global minima for N = 37 and 39. 

For LJ75 the global minimum was found in just one of 
the longer Monte Carlo runs, and again in the short run 
from the global minimum for LJ76. However, the latter 
minimum was only found in the short runs seeded from 
the LJ75 and LJ77 decahedra. Similarly, the LJ77 global 
minimum was only found in the short run seeded from 
the LJ76 decahedron. The decahedral global minimum 
for LJ75 was found in four out of 100 Monte Carlo runs 
of 5000 steps each, a frequency which fits in quite well 
with our results for LJ75, LJ76 and LJ77. A successful run 
requires an initial geometry which falls within the dec- 
ahedral catchment area; all the other runs produce the 
lowest icosahedral minimum after which the decahedron 
is never found. It would obviously be possible to locate 
global minima based upon decahedra more efficiently by 
biasing the starting configuration, but our intention was 
to analyze the performance of an unbiased algorithm in 
the present work. 

The pattern for LJ75-LJ77 is repeated for LJ102-LJ103. 
For LJ102 the decahedral global minimum was located in 
one of the longer Monte Carlo runs and in the short run 
seeded from the global minimum of LJ103. The decahe- 
dral minima for LJ103 and LJ104 were only found in short 
runs seeded from larger or smaller decahedra. The deca- 
hedral global minimum for LJ102 was found in three out 
of 100 Monte Carlo runs of 5000 steps each. 

The three new icosahedral global minima all have an 
atom missing from a vertex of the underlying Mackay 
icosahedron (Figure ||) . This is a possibility that Northby 
did not consider in his restricted search of the icosahedral 
configuration space. The new LJgg global minimum was 
found in three of the five longer Monte Carlo runs and 
in the short run seeded from LJ70. The new global min- 
imum for LJ78 was only found in the short run seeded 
from LJ79. The new minimum for LJ107 was found in 
one of the longer Monte Carlo runs and in the short run 
seeded from LJios- 

We also performed a few preliminary runs for LJ192 
and LJ201, sizes at which a complete Marks decahe- 
dron and a complete truncated octahedron occur, re- 
spectively. For LJ192 the Marks decahedron has energy 
-1175.697144. This structure was not found in 100 MC 
runs of 10000 steps each; instead the lowest minimum 
located had an energy of —1174.919801. For LJ 2 oi the 
truncated octahedron has energy —1232.731497. How- 



ever, we located a structure of energy —1236.124253 
which is based upon icosahedral packing. This minimum 
was found in three out of 50 MC runs of 10000 steps each. 
For these larger systems greater efficiency could proba- 
bly be achieved by varying the temperature and other 
parameters of the MC search. 

IV. CONCLUSIONS 

We have presented the results of a 'basin-hopping' or 
Monte Carlo-minimizationCJ approach to global opti- 
mization for atomic clusters bound by the Lennard- Jones 
potential containing up to 110 atoms. All the lowest 
known minima were located successfully, including the 
seven structures based upon fee or decahedral packing 
and three new global minima based upon icosahedra. Of 
the latter ten structures, only the smallest has been lo- 
cated before by an unbiased algorithm, to the best of our 
knowledge. 

The method is based upon a hypersurface deformation 
in which the potential energy surface (PES) is converted 
into a set of plateaus each corresponding to a basin of 
attraction of a local minimum on the original PES. This 
process removes all the transition state regions but does 
not affect the energies of the minima. On the original 
PES, most trajectories that approach the boundary be- 
tween two basins of attraction are reflected back due to 
the high potential energy; only if the trajectory is along 
a transition state valley does passage between basins be- 
come likely. In contrast, on the transformed PES it is 
feasible for the system to hop between basins at any 
point along the basin boundary which dramatically re- 
duces the time scale for interbasin motion. We speculate 
that the success of a previous genetic algorithm applied 
to the same clusters may be at least partly due to the 
fact that the same space is implicitly considered in that 
approach. E3 

The efficiency of the present approach could doubtless 
be improved by combining it with various other tech- 
niques. The most obvious short-cut would be to start 
not from initial random configurations but from seeds 
with either decahedral, icosahedral or fee morphologies. 
We have already checked that such biasing is indeed ef- 
fective, but our aim in the present paper was to gauge 
the performance of the unbiased algorithm. The temper- 
ature at which our Monte Carlo runs were conducted was 
also not optimized systematically. 

Finally, as we noted in the introduction, global op- 
timization for Lennard- Jones clusters at most sizes is 
a relatively easy task. A more stringent and general 
test is provided by Morse clusters which exhibit differ- 
ent structural-behaviour as a function of the range of 
the potentialoo At short range the task is particularly 
difficulL-because the PES is very xiigged — the number of 
minimaE£l and the barrier heightsEj increase as the range 
is decreased. 
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TABLE I. Global minima of LJat for N < 110. The references in which each minimum was first reported (to the best 
of our knowledge) are given, and t indicates the present wo rk. We intend to maintain an updated database of energies 
and coordinates for LJ and Morse clusters on our web site: [http://brian.ch.cam.ac.uk 



N 


Point group 


Energy /e 


Ref. N 


Point group 


Energy / e 


Ref. 


2 




-1.000000 


2 


57 


C s 


-288.342625 


10 




3 


D 3h 


-3.000000 


2 


| 58 


Czv 


-294.378148 


ic 




4 


T d 


-6.000000 


2 


3 59 


Civ 


-299.738070 


ic 




5 


D 3h 


-9.103852 


2 


3 60 


C s 


-305.875476 


ic 




6 


o h 


-12.712062 


2 


3 61 


Civ 


-312.008896 


ic 




7 


D 5h 


-16.505384 


2 


3 62 


C s 


-317.353901 


u 




8 


C s 


-19.821489 


1 


| 63 




-323.489734 


u 




9 


Civ 


-24.113360 


2 


3 64 


c s 


-329.620147 


u 




10 


C'iv 


-28.422532 


2 


65 


C-> 


-334.971532 


y 


'L_J 


11 


Civ 


-32.765970 


2 


66 


Ci 


-341.110599 


u 


12 




-37.967600 


2 


67 


C s 


-347.252007 


u 


13 


Ih 


-44.326801 


2, 


1 68 


Ci 


-353.394542 


iT 




14 




-47.845157 


2, 


| 69 


Csv 


-359.882566 






15 


Civ 


-52.322627 


2 


70 


C$v 


-366.892251 


10 




16 


Cs 


-56.815742 


2 


71 


Csv 


-373.349661 


ic 




17 


c 2 


-61.317995 


7 


72 


Cs 


-378.637253 


u 




18 




-66.530949 


2 


73 


Cs 


-384.789377 


17 




19 


D 5h 


-72.659782 


2 


74 


C s 


-390.908500 


i? 




20 


Civ 


-77.177043 


2 


75 


D 5h 


-397.492331 


u 




21 


Civ 


-81.684571 


2 


76 


Cs 


-402.894866 


~ 




22 


C s 


-86.809782 


m 77 


Civ 


-409.083517 


u 




23 


D 3h 


-92.844472 


s" 


78 


Cs 


-414.794401 






24 


Cs 


-97.348815 


9 


79 


Civ 


-421.810897 


E 




25 


C s 


-102.372663 


2 


80 


Cs 


-428.083564 


21 




26 


T d 


-108.315616 


2 


81 


Civ 


-434.343643 


iE 




27 


Civ 


-112.873584 


1( 


1 82 


Ci 


-440.550425 


iE 




28 


Cs 


-117.822402 


IT 


| 83 


Civ 


-446.924094 


21 




29 


D 3h 


-123.587371 




84 


Ci 


-452.657214 


£ 




30 


Civ 


-128.286571 


1C 


85 


Czv 


-459.055799 


iE 




31 


C a 


-133.586422 


k 


86 


Ci 


-465.384493 


IE 




32 


Civ 


-139.635524 


IC 


87 


C s 


-472.098165 


iE 




33 


Cs 


-144.842719 


1C 


88 


C s 


-479.032630 


iE 




34 


Civ 


-150.044528 


K 


89 


Czv 


-486.053911 


iE 




35 


Ci 


-155.756643 


K 


90 


Cs 


-492.433908 


iE 




36 


C s 


-161.825363 




91 


C s 


-498.811060 


iE 




37 


Ci 


-167.033672 


H 


92 


Czv 


-505.185309 


iE 




38 


o h 


-173.928427 


U 


,0 93 


Ci 


-510.877688 


iE 




39 


C§v 


-180.033185 


K 


94 


Ci 


-517.264131 


iE 




40 


Cs 


-185.249839 


H 


95 


Ci 


-523.640211 


iE 




41 


C s 


-190.536277 


ic 


96 


Ci 


-529.879146 


iE 




42 


Cs 


-196.277534 


K 


97 


Ci 


-536.681383 


iE 




43 


C s 


-202.364664 


K 


98 


C s 


-543.642957 


iE 




44 


Ci 


-207.688728 


K 


99 


Civ 


-550.666526 


iE 




45 


Ci 


-213.784862 


u 


100 


C s 


-557.039820 


iE 




46 


Civ 


-220.680330 


IT 


101 


C2v 


-563.411308 


>E 




47 


Ci 


-226.012256 


IT 


102 


C2v 


-569.363652 


>E 




48 


C s 


-232.199529 


IT 


103 


Cs 


-575.766131 


iE 




49 


c 3v 


-239.091864 


IT 


104 


C2v 


-582.086642 


iE 




50 


Cs 


-244.549926 


IT 


105 


Ci 


-588.266501 


iE 




51 


C2v 


-251.253964 


IT 


106 


Ci 


-595.061072 






52 


c 3v 


-258.229991 


IT 


107 


C s 


-602.007110 


T 


53 


C2v 


-265.203016 


IT 


108 


Cs 


-609.033011 


Ifj 


54 


C$ v 


-272.208631 


IT 


109 


c\ 


-615.411166 


ni 


55 


Ih 


-279.248470 


| 


110 


C a 


-621.788224 


itj 


56 


c 3v 


-283.643105 


k 











G 



TABLE II. Lowest energy icosahedral minima at sizes 
with non-icosahedral global minima. 



IN 


Point group 


Energy / e 


Ref. 


38 


c$ v 


-173.252378 




16 


75 


Ci 


-396.282249 




n 


76 


Ci 


-402.384580 




12 


77 


Ci 


-408.518265 




12 


102 


C s 


-569.277721 




io 


103 


Ci 


-575.658879 






104 


C s 


-582.038429 
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